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ABSTRACT 



For a large class of dark energy (DE) models, for which the effective gravitational 
constant is a constant and there is no direct exchange of energy between DE and dark 
matter (DM), knowledge of the expansion history suffices to reconstruct the growth 
factor of linearized density perturbations in the non-relativistic matter component 
on scales much smaller than the Hubble distance. In this paper we develop a non- 
parametric method for extracting information about the perturbative growth factor 
from data pertaining to the luminosity or angular size distances. A comparison of the 
reconstructed density contrast with observations of large scale structure and gravita- 
tional lensing can help distinguish DE models such as the cosmological constant and 
quintessence from models based on modified gravity theories as well as models in which 
DE and DM are either unified, or interact directly. We show that for current SNe data, 
the linear growth factor at z = 0.3 can be constrained to 5%, and the linear growth rate 
to 6%. With future SNe data, such as expected from the JDEM mission, we may be 
able to constrain the growth factor to 2 — 3% and the growth rate to 3 — 4% at z = 0.3 
with this unbiased, model-independent reconstruction method. For future BAO data 
which would deliver measurements of both the angular diameter distance and Hubble 
parameter, it should be possible to constrain the growth factor at z = 2.5 to 9%. 
These constraints grow tighter with the errors on the datasets. With a large quantity 
of data expected in the next few years, this method can emerge as a competitive tool 
for distinguishing between different models of dark energy. 

Subject headings: cosmology: cosmological parameters — cosmology: distance scale 
- cosmology: theory 
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Introduction 



Over the last decade, observations of Type 



Perlmutter et al. 1999 




Tonrv et al. 2003 


Riess et al. 2007: 


Wood-Vasev et al. 2007 



a supernovae have shown that the ex 



Perlmutter et al. 1998: 



Riess et al. 2005 



Riess et al. 1998 



Astier et al. 2005 



Kowalski 20081 ). This remarkable discovery has 



led cosmologists to hypothesize the presence of dark energy (DE), a negative pressure energy 
component which dominates the energy content of the universe at present. Many theories 
have been propounded to explain this phenomenon, the simplest of which is the cosmological 
constant A, with a constant energy density and the equation of state w = —1. Although A 
appears to explain all current observations satisfactorily, to do so its value must necessarily 
be very small A/8ttG ~ 10~ 47 GeV 4 . So, it represents a new small constant of nature in 
addition to those known from elementary particle physics, many of them being very small 
if expressed in the Planck units. However, since it is not known at present how to derive 
A from these small constants and it is also unclear if DE is in fact time independent, other 



Sahni & Starobinskv 2000 




Carroll 2001 




Pee : 


les & Ratra 2003 




Padmanabhan 2003 


Copeland et al. 2006: 


Noiiri & Odinstov 2007: 


Sahni & Starobinskv 2006 


). These are 



either on the introduction of new physical fields (quintessence models, Chaplygin gas, etc.), 
or on modifying the laws of gravity and therefore the geometry of the universe (scalar-tensor 
gravity, f(R) gravity, higher dimensional 'Braneworld' models e.t.c. ). The plethora of 
competing dark energy models has led to the development of parametric and non-parametric 
methods as a means of obtaining model in dependent information about the nature of 



dark energy directly from observations (see 



Corasaniti et al. 2003 



Starobinskv 1998: 



Huterer fc Starkman 2003 



Huterer fc Starkman 1999 



Alam et al. 2004 



Saini et al. 2004 



Jassal et al. 2005 



Wang fc Mukherjee 2006 



Lazkoz et al. 2006 



Alam et al 2007 



Sarkar et al. 20081 : ISahni et al. 20081 : ISahni fc Starobinskv 20061 . and references therein) . 
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The next decade will see the emergence of many new cosmological probes. A large 
number of these are likely to make important contributions to the field of dark energy. The 
Sloan Digital Sky Survey began its stage III observations in 2008, and its Baryon Oscillation 
Spectroscopic Survey (BOSS) is expected to map the spatial distribution of luminous 
galaxies and quasars and detect the characteristic sca le imprinted by baryon acoustic 



oscillations in the early universe (ISDSS collaboration! ) . The Joint Dark Energy Mission 



(JDEM) is expected to discover a large number of sup ernovae, and also prov ide important 
data on weak-lensing and baryon acoustic oscillations (IJDEM collaboration! ) . The Square 
Kilometer Array (SKA) will map out over a billion galaxies to redshift of about 1.5, and is 
expected to determine the po wer spectrum of dark matter fluctuations as well as its growth 



as a function of cosmic epoch (IBlake 20041 ). Important clues to the growth of structure will 



also come from current and future weak lensing surveys (CFHTLS, PES, JDEM, EUC LID 



SKA, LSST), galaxy redshift-space distortions (IGuzzo 2008 



Song 2008 



Percival 2008T) as 



Rapetti 2008h . 



well as galaxy cluster mass functions at different redshifts z ( jVikhlinin 2008 
With the wealth of data expected to arrive over the next several years, it is important 
to explore different methods of analyzing these datasets in order to extract the optimum 
amount of information from them. In this paper we explore the possibility of reconstructing 
the linearized growth rate of density perturbations in the non-relativistic matter component, 
S(z), taken at some fixed comoving scale much less than the Hubble distance, from datasets 
which have traditionally been used to explore only the smooth background universe, 
e.g. luminosity distance and angular diameter distance data. 

In the case of physical DE, the effective gravitational constant appearing in the 
equation for linear density perturbations in the matter component coincides with the 
Newton gravitational constant G measured in the laboratory and using Solar system 
tests. If, additionally, there is no direct non-gravitational interaction between DE 



and DM in the physical reference frame, so that the DE energy-momentum tensor is 
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covariantly conserved, the density contrast reconstructed in this manner should match 
that determined directly from observations of large scale structure. In this case the 
methods developed in this paper will provide an important consistency check on DE 
models such as A and Quintessence. On the other hand, geometrical models of DE 
(Braneworlds, scalar-tensor gravity, etc.) usually predict a different growth rate for 5(z) 
from that in general relativity (GR). Models where DE has a direct non-gravitational 
interaction with DM, or where DE and DM are unified, have a similar property 
even in the framework of GR. In this reconstruction of the linearized density 

contrast from observations of stan dard candles /rulers will not match with 5 determinec 



directly from large scale structure ( 



Knox et al. 2006 



Lue et al. 200 



Chiba fe Takahashi 2007 



Bertschinger 20061: 



Huterer Sz Linder 2007 



Ishak et al. 2006 



Wang et al. 20071 : 



Mortonson et al. 2008 



Polarski fc Gannouji 2008 



Song fc Koyama 20081 ) . Currently 



reconstructed va 



Verde et al. 2002 



Hawkins et al. 20Q[ 



future missions like Euclid ( jCimatti et al. 20081 ) are expected to constrain the growth rate 



u es of the growth rat e from galaxy red shift distortions (IGuzzo 2008 
"J 



Ross et al. 20071 ) are not very constrictive, but 



tightly. Therefore comparing the results from future supernova data, using the methods 
described in this paper, to those from future large scale structure data will help address 
important issues concerning the nature of gravity and dark energy. 

This paper is organized as follows. In section [2], we describe the reconstruction 
technique and the data used to test this method. Section [3] shows the results and examines 
the dependence of the method on various factors such as the redshift distribution of the 
data and information on other cosmological parameters. The conclusions are presented in 
section HI 
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2. Methodology 



In the longitudinal (quasi-Newtonian) gauge, the perturbed, spatially flat, Friedman- 
Robert son- Walker (FRW) metric is defined by the line element 



1 + 2<p)dt 2 + (1 - 2ij)a 2 (t)dx 2 



where = if) in GR if matter is free of anisotropic stresses (we assume DM to be cold and 
neglect small effects from the neutrino component which produces (f> ^ if>). The Newtonian 
potential <fi and the non-relativistic matter density contrast 

. Pm(x,t) ~ p(t) 



P(t) 



(2) 



are linked via the linearized Poisson equation 



k 2 cf) = —AnGa 2 p m 5 ri 



If the DE energy-momentum tensor is covariantly conserved, then 



(3) 



case it is straightforward to show (see, e.g. 



Wang & Steinhardt 1998 



p m oc a 3 . In thi s 



Starobinskv 19981 ) 



that on scales much smaller than the effective Jeans scale for DE, Aj ~ CgH" 1 , where 
H(z) = a/a is the Hubble parameter and c s is the effective DE sound velocity (c s = 1 
for standard quintessence), linearized matter density perturbations in a FRW universe 
containing DE with an arbitrary effective equation of state w(t) = Pde/pde satisfy the 
same equ ation as in the case of a standard FRW model driven by dust and a cosmological 



constant (IPeebles 19801 ) : 



2H5 m - AnGp m 5 m = 



(4) 



we igno re the subscrip 



solution (IPeebles 1980 



; in S m in the ensuin g discussion). However, the generic textbook 



Sahni fc Coles 19951 ) 



5 oc H(z) 



•-0 



1 + 31 



dzi 



(5) 
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is not applicable now, apart from the following cases: dust-like matter, a non-zero spatial 
curvature (and/or a tangled network of cosmic strings) and a cosmological constant, for 
which H 2 (z) = C\ + 6*2(1 + z) 2 + 6*3(1 + z) 3 . The same refers to the other well-known 
expression valid only for dust and a cosmological constant: 



8 oc a(t) - H(t) j a{t x )dti, (6) 



see e.g. ( jKofman fc Starobinsky 1985 



Bertschinger 20061 ) J 1 ! Thus, for an arbitrary physical 



DE, Eq (jlj) has to be solved numerically. Since there are no terms depending on the 
perturbation wave vector k in it, 8(z)/5(0) will be fc-independent, too. We will also suppose 
that c s is not too small, so that k ^> a/Xj for all scales of interest, in particular, c s ^> 0.01 
if we consider scales up to 100(1 + z)^ 1 Mpc. 
The dimensionless physical distance 

where to is the present moment, plays a key role in measurements of the background universe 



x The expression is, in fact, the first term in the long-wave (super-Hubble) expansion 
of the adiabatic mode of a comoving density perturbation if the perturbed matter pressure 
tensor is proportional to the unit one and the spatial curvature may be neglected. Then 
(j) — ip — —((x) (l — ^ J 4 * adt^j, where the gauge-invariant c urvature perturbation C does 



not depend on tim e for the growing adiabatic mode, see e.g. (IPolarski fc Starobinsky 1992 



1 

Bertschinger 20061 ). The quantity ti is free and may be chosen to coincide with the moment of 



the first Hubble radius crossing during inflation (another choice would correspond to adding 
a decaying adiabatic mode with an arbitrary amplitude). Then, using Eq. ([3]) and the fact 
that p m oc a" 3 , the formula (|6]) follows. Since DE is practically unclustered at sub- Hubble 
scales, it is tempting to try to use this formula for Aj A <ti H~ l , too. However, as pointed 
above, this works only if DE is a cosmological constant. 
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using standard rulers and candles. E is proportional to the conformal time measured from 
the present to a moment in the past. It is related to the luminosity distance, cIl, via 



H d L (z) 
1 + z 

H d L (z) 
1 + z 

H d L (z) 
l + z 



1 



E(z) , 
1 



sm{^/\Q^\E(z)} 



smh{^n K E(z)} 



n K < o 

tt K = 
tt K > 



(9) 
(10) 



K 



where fix = 1 — ^totai- The following relationship between the luminosity distance d^ and 
the angular size distance d^ holds in a metric theory of gravity: — (1 + z) 2 dj±. Rewriting 
Eq (j3J) in terms of Eq j7]) and using the fact that p m oc (1 + z) 3 , we obtain: 



5' 



-Qo m 5 



11) 



l + z(E) 

where the prime denotes a derivative with respect to E. It is straightforward to transform 
a (1111) into the following set of integral equations for 5(E) and its first derivative 



Sahni fc Starobinskv 20061 ) 



5(E) = l + 5' [l+z(E 1 )]dE 1 + -Q 0m [l+z(Ex)] 
Jo * Jo 

5'(E) = 5 [l + z(E)] + ~n 0m [l + z(E)} [ 5(E 1 )dE 1 , 

* Jo 



5(E 2 )dE 2 dEx (12) 



(13) 



where 5 is normalized to 5q = 5(z = 0) = 1. Note the remarkabl e fact that, in 



contrast to formulas used in the re construction of H (z) from d L (z) (jStarobinsky 1998 



Huterer Starkman 19991 ) or 5(z) (jStarobinsky 19981 ) which require taking a derivative 
of observational data with respect to the redshift, this formula contains integrations of 
observational data only, which is a sound operation for noisy data. 

By solving the above equations we can calculate the linear growth factor 



g(z) = (l+z)5(z) 



(14) 
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which represents the ratio of 5(z) in the presence of dark energy to that in SCDM without 
a cosmological constant. Another quantity of interest is the growth rate 

m = £i = -111 m . (15) 

a ma H(z) o(E) 

To solve Eq ffl2l) we start with initial guess values for 5(E) and 5'(E) and iteratively 
solve for 5(E), calculating 5'(E) in the successive iterations as the difference between 
adjacent values of 5(E), i.e. 5[ = A5i/AEi. For calculating f(z) we require to estimate 
H(z) as well. We obtain this quantity by differentiating the noisy data E(z) using a finite 
differencing method . This naturally amplifies the noise in the final results, so we expect the 
results for f(z) to be slightly noisier. However, typically the difference in f(z) between two 
models of dark energy is greater than the difference in g(z), so despite the greater noise, 
we expect f(z) to be useful for discriminating dark energy models. This method does not 
require prior knowledge of the parameter 5' Q , is robust to changes in the initial guess values 
and gives exact results for g(E) and f(E) for noiseless data. For data with errors, naturally 
the result is noisier, however, as we will show in the succeeding sections, we will be able to 
put reasonable constraints on g and / using this method. 

Data noise can also be decreased using smooth ing techniques. In wh at follows we shall 



use the lognormal smoothing scheme proposed in (jShafieloo et al. 20061 ) which has been 
shown to be reasonably unbiased and efficient. It const ructs a smooth quantity, E s , from a 
noisy one, E(zi), via the ansatz (jShafieloo et al. 20061 ) 



l n 2i±2A _ /-ln 21+ ^ 



E>{z) = X>(*)exp ( /5>P ) , (16) 



where A is the smoothing scale (see also 



Shafieloo 20071 ). We take A ~ 1/N where N is the 



total number of observations. Choosing this small value of A leaves the results unbiased. 

From the manner in which Eqs fl X lti X 3 j) have been obtained, reasons as to why the 
linearized growth function 5 bs(z) determined from actual observations of large scale 
structure may differ from 5(z) reconstructed using our method follow immediately: 
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1. p m is not proportional to (1 + z) 3 . This happens even in GR if the energy- momentum 
tensor of physical DE is not covariantly conserved separately, either due to the existen ce 
of a direct non-gravita t ional interaction between DE and DM (see 



Billvard fc Colev 2 000 



Zimdahl k Pavon 2001 



Amendola 1999 



Caldera-Cabral et al. 2008 



, and 



references therein), or because DE and DM constitute some unique entity as oc curs 
in unified DM-DE models such as the Chaplygin gas (IKamenshchik et al. 20011 ) and 
its generalizations. In such models DE is partially clustered with DM on small scales, 
and this can result in the appearance of significant fc-dependent terms in the equation 
for S, so that the growth factor g(z) becomes fc-dependent. In particular, the latter 
effect is especially crucia l for the generalized Chaplygin gas model, (see the recent 



paper 



Gorini et al. 20081 . and references to previous papers therein). 



2. GR is modified, DE is geometrical. Then G in Eq (J4]) becomes some effective quantity 
G e ff which may be both time and scale dependent. A noticeable value of — if) may 
also arise even in the absence of free-streaming particles. However on small scales this 
value is strongly restricted by Solar system tests of gravity (which don't suggest any 
such effect). Geometrical models of DE, which include Braneworld models and models 
using scalar-tensor and f(R) gravity, have more degrees of freedom than GR, so it is 
natural that in this case, the linearized perturbation equation for S shows a departure 
from the Newtonian form, Eq flU). For in stance in extra dimensional scenario's 



flDvali et al. 200C 



Sahni &: Shtanov 20031 ). the presence of the fifth dimension 



(the bulk) can influence the behaviour of perturbations residing on the brane 



flLue et al. 2004 



Koyama fc Maartens 2006 



Sawicki et al. 2007 



Shtanov et al. 20071 ) 



making t hem significantly fc-d ependent even on scales much smaller than the Hubble 
distance l 1 



gravity 



review 



Shtanov et al. 20071). The same effect arises in viable DE m odels in f(R) 
Hu fc Sawicki 20071 ; IStarobinsky 20071 ; iTsujikawa et al. 20081 ) (see the recent 



Sotirou &: Faraoni 2008 



for numerous papers on f(R) gravity as a whole). 
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On the other hand, in some cases G e fj and g(z) may remain scale-independent even 
on small scales, though they acquire a non-trivial, non-GR, time dependence. This 
occurs, for instance, in the Dvali-Gabadadze-Porrati (DGP) extra-dimensional model 
considered below as well as in sc alar-tensor DE models with a large current value of 



the Brans-Dicke parameter ubd ( jBoisseau et al. 2000 



Gannouji fc Polarski 20081 ) 



Thus, a comparison of the observed and reconstructed density contrast could help 
shed light on the nature o f dark energy. While it is encou raging that future observations 



( IBlake 2004 



Guzzo 20081 : ISong 20081 : ICimatti et al. 20081 1 of large scale structure may 



make possible the determination of S b s (z), in this paper we focus on reconstructing 5 using 
observations of high redshift type la supernovae and baryon acoustic oscillations. 



2.1. Data used 

The method outlined in the previous section would be applicable to any observation 
which contains a measurement of E(z), e.g. measurements of luminosity distance or angular 
diameter distance. We shall use real data and mock data based on simulations of supernova 
type la data and the angular diameter distance from baryon acoustic oscillations, to test 
this method. 
Supernova Data : 

The lightcurves of Type la supernovae show them to be "calibrated candles" , therefore 
they are of enormous significance in cosmology today. The luminosity distance of Type la 
SNe provide us with a direct measurement of the acceleration of the universe, thus leading 
to constraints on the dark energy parameters. SNe data is in the form {mg, z, <J mB , a z }, 
where the magnitude rne is related to cIl(z) as 

m B = 5log w [H d L (z)] +M, (17) 

M. being a noise parameter usually marginalized over. 
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It should be noted that supernova data, on its own, is unable to break the degeneracy 
between dark energy and spatial curvature. The CMB, on the other hand, places stringent 
constraints on Qx and strongly suggests that the universe is spatially flat, in agreement 
with predictions made by the inflationary scenario. In this paper we shall work under the 
assumption that fix = and use to relate di — > E(z), with the latter playing the key 
role in our reconstruction exercise (lT2"j) . 

Currently there are around 30 published SNe with the furthest observed one at a 
redshift of z = 1.7 (IKowalski 20081 ). and average error o f a mp ~ 0.15. Future space-based 



projects such as the Joint Dark Energy Mission (JDEM) (IJDEM collaboration! ) are expected 



to observe about 2000 SNe with errors of a 



m B 



0.07. To date, SNe are the most direct 



evidence for dark energy, and in this paper we shall primarily use SNe data to constrain the 
growth parameters for different dark energy models. 
BAO data 

At present, baryon acoustic oscillations are believed to be the method least plagued 
by sys tematic uncertainties, therefore the detection of the first baryon acoustic oscillation 



scale (jEisenstein et al. 20051 ) has led to the speculation that BAO may in future become 
a potent discriminator for dark energy. Standing sound waves that propagate in the 
opaque early universe imprint a characteristic scale in the clustering of matter, providing 
a "standard ruler". Since the sound horizon is tightly constrained by cosmic microwave 
background (CMB) observations, measuring the angle subtended by this scale determines 
a distance to that redshift and constrains the expansion rate. The radial and transverse 
scales give measurements of [r s H(z)]/c and r a /[(l + z)cIa(z)} respectively, where r s is the 
sound horizon obtained from CMB. These quantities are correlated, and the present BAO 
data is not se nsitive enough to measure both quantities independently (see however the 



recent papers 



Benitez et al. 2008 



to give independent measurements of cIa(z) and H(z) (jSeo fc Eistenstein 20031 ). Future 



Gaztanaga et al. 2008T). but future surveys a re expected 
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BAO surveys such as BOSS (jSDSS collaboration! ) should therefore place tighter constraints 
on dark energy parameters. 

3. Results 

We first use Supernova data to reconstruct the growth parameters. We simulate data 
according to two theoretical models : 

• Model 1 : A cosmological constant model with w = — l,f^om — 0.27, H = 
72 km/s/Mpc. 

• Model 2: A variable dark energy model with the equation of state given by 

vj z 

w(z) = w + —— , w = -0.9, w a = 0.3 , (18) 
i. ~\~ z 

and with the same values of f2o m , H as Model 1. Note that the Models 1 and 2 provi de 



excellent agreement with the current CMB+BAO+SNe data (IKomatsu et al. 20081 ). 
The Model 2 has w > — 1 everywhere, so it can be realized by quintessence with some 
potential for which Eq (jl]) is valid. 

Two different data distributions are used, set A resembles the quality of data available 
at present, and set B is modeled on expected future surveys. 

• Set A : ~ 300 SNe, with the redshift distribution and errors of the Union dataset 



( iKowalski 20081 ). For this dataset, on average, a mB ~ 0.15, but a few SNe have very 
high errors of the order of unity. Since the method of integration would not work 
very well for very noisy data, and a single datapoint with large noise would affect the 
results of all datapoints after it, we restrict the analysis to SNe with a mB < 0.7. By 
rejecting only 10 datapoints with this criterion, we enhance the results by a significant 
amount. 
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Set B : ~ 2000 SNe, with the redshift di stribution and errors [p mB ~ 0.07) expected 



from future surveys such as the JDEM (jAldering et al. 20041 ). The data covers a 



redshift range of z — 0.1 — 1.7 with a larger concentration of supernovae in the 
midrange redshifts (z = 0.4 — 1.1). The errors considered here are statistical only, we 
do not consider systematic errors, which are expected to be better controlled in the 
future with larger datasets. 

For both cases, we marginalize over Q Qm = 0.27 ± 0.03. Supernova data is unable to break 
the degeneracy between dark energy and curvature of the universe. In order to measure the 
growth parameters, we therefore consider only a flat universe, which is the preferred model 
from current CMB observations. 

Fig [T] shows the results for the linear growth factor g(z) for both datasets and for the 
two different cosmological models. We see that for both models, set A results in rather noisy 
reconstruction (left panel), since the errors on the SNe are quite high. This is especially 
true at high redshifts (z > 0.7) where the sparse sampling affects the integral reconstruction 
scheme adversely. At z = 0.3, g(z) is constrained to ~ 5%. For JDEM-like data (set B) 
however, g(z) is reconstructed more accurately, and has low errors at low redshifts (right 
panel). At z — 0.3, g(z) is constrained accurately to ~ 2% for both models for set B, while 
at z — 1, g(z) is constrained to ~ 4%. 

Fig [2] shows the reconstruction of the growth rate f(z). As before, the results for set A 
are poor, with f(z) constrained to ~ 6% at z — 0.3. The results for Set B are reasonable, 
however, the errors are slightly larger in this case than for g(z), since there is an additional 
error from the calculation of H(z) from E(z). At z — 0.3, f(z) is constrained accurately to 
~ 3% for both models for set B, while at z — 1, f(z) is constrained to ~ 8%. We also note 
that the quantity f(z) has slightly greater discriminatory power than g(z), since typically 
the growth factor shows rather less variation between different dark energy models as 
compared to the growth rate at any given redshift. Therefore, even though f(z) is slightly 
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noisier, for set B, Model 1 and Model 2 can be discriminated at la using f(z). 

If the data is first smoothed with the smoothing scheme, Eq (|T6|) . the results improve, 
especially for set A which has much noisier data, as seen in figure [31 The results for 
f(z) improve markedly for both datasets. This is because an additional quantity H(z) is 
required for obtaining f(z), and a smoother E(z) leads to a much more accurate estimation 
of H(z). Errors on g(z) and f(z) are ~ 1% and ~ 1.5% respectively at z — 0.3, and ~ 3% 
and ~ 6% respectively at z — 1 for Model 1 with JDEM like data. Model 2 gives similar 
constraints. The results for the growth parameters are summarized in Tabled] for Model 1, 
and in Table [2] for Model 2. We see that this method obtains quite reasonable constraints 
on the growth parameters at low redshifts for the set B, therefore it can be used successfully 
to constrain growth parameters from future SNe data. It should be noted that, for future 
SNe data to accurately constrain the growth parameters, it is important to keep the SNe 
systematics under control (a sys < 0.05). A systematic error of a sys = 0.1 (as on the current 
data) would weaken all constraints significantly. 

3.1. Dependence on nature of data 

We now check how the results change if the redshift distribution or error distribution 
is changed. To study the dependency on the number of SNe, we use three redshift 
distributions- (a) set A (~ 300 SNe) with double the number of supernovae at low (z < 0.3) 
and high(z > 0.7) redshifts, (b) set A with double the SNe at mid-range (0.3 < z < 0.7) 
redshifts, and (c) a distribution with the JDEM (set B) redshift distribution (~ 2000 SNe) 
with errors of the order of the Union (set A) SNe. The results for Model 1 are shown in 
figure HI We see that doubling the number of SNe in a particular redshift bin changes 
the results very slightly. This is to be expected because when integrating noisy data, 
having a larger number of points with the same amount of noise does not improve results 
significantly. Increasing the total number of SNe by a significant amount (nearly seven 



-16- 



times, as in right panel) does improve the scatter, but the results still do not compare with 
those of set B (fig [21 top right panel) which has the same number of supernovae but smaller 
errors. 

We now study the effect of the errors. Once again we study three distributions — (a) 
set A with the errors halved for z < 0.3 and z > 0.7 redshift bins, (b) set A with errors 
halved in the 0.3 < z < 0.7 redshift bin, and (c) set A with errors replaced by JDEM-like 
errors on all SNe. The results for Model 1 are shown in figure [5j We see that in this case, 
decreasing the errors at low redshift or high redshift changes the results very slightly. This 
is because there are very few points at low redshift so they do not affect the integration 
process strongly, while the high redshift points cannot affect the low redshift points. The 
results in the redshift range 0.3 < z < 0.7 become better if the mid-range SNe have lower 
errors. As we see in the right panel of fig El decreasing the errors to JDEM errors gives 
results almost identical to the results for Set B (fig El top right panel) , even though the 
number of points is much less for set A. Thus we find that this method would work quite 
well even for a reasonable number of supernovae (of the order of a few hundred) provided 
the errors were tightly constrained. 

Since the high errors of set A make it unsuitable for this reconstruction approach, in 
the next sections we will use the set B to study the robustness of the results to various 
other factors. 



3.2. Growth rate from w(z) 



We may also calculate the growth r ate / from the supernova data via the equation of 



state using the following approximation (jWang fc Steinhardt 19981 ) : 



/(*) 
j(z) 



flom(l + Z) 

H 2 (z) 



5- 



+ 



3 (l-iu)(l-§to) 



l—w 



125 



;i - iwf 



{i-a m {z)) + o\{i-n m {z))f 



(19) 
(20) 
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where the equation of state w(z) may be calculated using a likelihood parameter estimation 
from the luminosity distance. This approximation works quite well for a large number 
of physical dark ene rgy models with a constant or slowly c hanging w including LCDM, 



for which 7 ~ 0.55 (IWang fc Steinhardt 1998 



flChevallier fc Polarski 2001 



Linder 200.T 1 : 



Linder 20051 ). We use the familiar CPL fit 



w(z) 
H 2 (z) 



w a z 
w + -— , 
1 + z 



Ht[n 0m (l + z) 3 + (l-tt 0m )(l + z) 



3(l+u) +w (1 )„3w (1 (l/(l+^)-l) 



(21) 
(22) 



A likelihood parameter estimation is expected to lead to smaller errors, but the drawback 
of this method is that the result may be biased due to the parameterization. Also the 
errors on w(z) would propagate extremely non-linearly to / and therefore the result for 
f(z) would be much less trustworthy. 

Figure [6] shows the reconstructed f(z) for Model 1 and 2 for set B. As expected, the 
errors are lower that those for our reconstruction method. However, it is also noteworthy 
that the resulting confidence levels are not symmetric around the true value, in fact at 
higher redshifts, the true model appears to be on the verge of b eing ruled out ! These results 



are commensurate with those found in (IMortonson et al. 20081 ). where reconstruction of the 



growth parameters through w leads to biases in the growth parameter results even though 
w is recovered accurately. This is due to the fact that errors propagate non-linearly from 
w to f(z). We therefore conclude that, when reconstructing the growth parameters from 
supernova data, it is better to reconstruct the quantities directly, rather than reconstructing 
them indirectly from the energy density or equation of state. 



3.3. Dependence on Oom 

Supernova data do not simultaneously constrain information on fl 0m and dark energy 
parameters. To reconstruct dark energy parameters, it is necessary to place constraints on 
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fiom from other observations. In the calculations so far, we have marginalized over the true 
fiducial value for f2 0m . However, since there is considerable uncertainty as to the real value 
of the matter density, we check how using incorrect values of Qo m may bias our analysis. 
(It is well known that an incorrect value of Qom can significantly bias the results for DE 



flShafieloo et al. 2006 



Sahni et al. 20081 ).) The fiducial universe for model 1 contains 
n 0m — 0.27. We now choose a different, incorrect, value of n m = 0.3 for marginalization 
and proceed to analyze the data using both the integral reconstruction method and the 
likelihood parameter estimation of w outlined in the previous section. The results are 
shown in figure H We see that choosing a higher value of n 0m gives biased results in both 
methods, but interestingly enough, the biases are in opposite directions ! In case of the 
integral reconstruction method, a higher value of f2om leads to a lower value of f(z) at high 
redshifts, whereas for the w parameterization, a higher value of f2om leads to a higher value 

of/ (4 

These results may be understood as follows. For the reconstruction from w, we see from 
eq (Tl9|) that f(z) changes primarily due to the change in the matter density fW(l + z ) 3 , 
since the value of 7 does not vary very strongly with w, and H 2 (z) is constrained by the 
data. Choosing a higher value of f2o m would result in the choice of a different w(z) which 
would lead to nearly the same H(z) as that for a lower value of f2om, and 7 would also not 
change by much. However, the quantity n 0m (l + z) 3 would increase proportionate to fi 0m . 
Therefore a higher value of Oom would simply result in a higher value of f{z). In the case 
of the integral reconstruction however, we see from eq ( fl2i) that both 5 and 5' depend on 
O 0m - In 8 the leading term is unity and the other two terms containing S' and fl 0m are at 
about an order of magnitude smaller. In 5' the two terms containing 5' and Qo m are of 
the same order and opposite sign. The Qo m term contributes by making 5' less negative. 
Therefore increasing Q 0m increases S slightly and decreases the absolute value of 5' by a 
larger amount, so that the ratio between 5 and 5' becomes a smaller negative quantity. 
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Since f(z) is essentially this ratio, this means that f(z) also decreases with increasing Qo m . 
Therefore, choosing a wrong value of f2 0m causes the two different methods of reconstruction 
to be biased in opposite directions. This leads to the interesting conclusion that, provided 
other systematics are under control, comparing the integral reconstruction method with the 
standard likelihood estimation would give us a valuable consistency check on the accuracy 
of the prior chosen for O 0m . 



3.4. Reconstruction for a toy modified gravity model 



An influentia l Braneworld model was suggested by Dvali-Gabadadze-Porrati 



( IDvali et al. 20001 ). The expansion history for this model is given by 

1 — ^Om 



Hz) 



Hn 



1-0, 



Om 



(23) 



For physical models of da rk energy, the growth rate is well approximated by ea (fT9l) . fo r 



instance -7 ~ 0.55 for LCDM flWang fc Steinhardt 1998 



Thomas et al. 2008 



Linder 2005 



Peacock et al. 2001 



Acquaviva et al. 20081 ). This equation is not valid however if the 



observed acceleration originates from a modification of the equations of general theory 
of relativity; in the DGP Braneworld theory, the growth rate is approximated by 



(ILue et al. 2004 ) 



f(z) ~ VL m (z 



,0.68 



(24) 



This is the growth rate which would be measured through galaxy redshift distortions or 
weak gravitational lensing, whereas any analysis from the expansion history would obtain a 
growth rate commensurate with eq ( fT9l) . 

Therefore, if the growth rate for this model is reconstructed using on the one hand, 
supernova data, and on the other, galaxy redshift distortions, we expect the results to be 
different. We reconstruct the growth rate using the JDEM-like SNe distribution for this 
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modified gravity model by substituting Eq fl23|) into the integral reconstruction method 
described by Eq f ll2p . The result is shown in figure [HJ For comparison, we have also plotted 
the expected o bservational constrai nts from galaxy redshift distortions for the future 



Euclid mission (ICimatti et al. 20081 ). We see that the two results are strongly discrepant, 
especially at low redshifts. If the origin of dark energy were indeed geometrical in nature, 
comparisons of this sort would provide crucial evidence for it. 

Despite the popularity of the DGP model, it is currently facing several difficulties both 



of an observational and theoretical nature: Tension between t 



lis model and observational 



data sets has been pointed out in (jFairbairn fe Goobar 2006 



Alam fc Sahni 2002 



gravity ( 



Alam fc Sahni 2006 



Maartens fc Maioretto 20061) and the presence of a ghost in DGP 



Charmousis et al. 20061 ; 



Gregory et al. 2007 



Deffavet et al. 2006 



Kovama 20071 ) 



may be even more problematic. Consequently our purpose in the present section has been to 
treat DGP cosmology as a toy model, used to demonstrate the utility of the reconstruction 
approach developed i n this paper. (Note however the existence of other braneworld models 



which are gh ost free ( 



Sahni fc Shtanov 2 003 
T 



Shtanov et al. 20091 1 and agree well with 



observations (lAlam fc Sahni 20061 ).) 



3.5. Current SNe Data 

In figure [9] we show the reconst ructed growth p arameters for the currently available 



supernova data- the Union dataset (jKowalski 20081 ) . The results are marginal ized over 



Oom = 0.26 ± 0.03, the currently accepted value of f2om ( IDunkley et al. 20081 ). The 
nuisance parameter M. which contains information on Hq is also marginalized over. For the 
non-smoothed method, since errors are quite large, it is difficult to put any constraints on 
the growth parameters. If the smoothing scheme is used, f(z) may be constrained to ~ 6% 
at z — 0.3. At this redshift, the growth factor g(z) would be constrained to ~ 5%. The 
reconstructed f(z) is commensurate with the cosmological constant model as well as Model 2 
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(variable w, eq (ITS1) ) used in this paper. We also show the three curren t observations of H z) 



from galaxy redshift-space distortions (IGuzzo 2008 



Verde et al. 2002 



Ross et al. 20071 ). 



The error bars on these observations are at present quite large, but it is expected that future 
data in this field would be comparable with our results from supernovae, thus we would be 
able to discern physical and geometrical DE using these different techniques (as shown in 
section I3"^E1) . Table [3] shows the la limits on the growth parameters for the reconstruction. 



3.6. Data expected from future BAO experiments 

We now check the method with BAO data. The SDSS baryon acoustic oscillation 
survey of BOSS is expected to measure the baryon acoustic oscillation power spectrum 
very accurately. The expected accuracy on the angular diameter distance (1a is of the order 
of 1.0% at z = 0.35, 1.1% at z = 0.6 , and 1.5% at z = 2 .5, with errors on H(z) of 1.8%, 



1.7% and 1.5% at the same redshifts (ISDSS collaboration! ) . We populate a redshift range of 



z = 0.2 — 2.5 with 20 datapoints with errors based on these numbers and use this dataset 
to reconstruct the growth parameters. Since there are only 20 points in the dataset, and 
not many at very low redshifts, the integration is not very accurate even though the errors 
on d,A and H are small. We find that for this dataset g(z) and f(z) are both constrained 
to ~ 9% at z — 2.5 (see Tables [TJ and [21 bottom row). Although these errors appear to 
be large compared to those from the SNe data, for a high redshift of z — 2.5, these errors 
are actually commensurate to the errors from SNe. The advantage of using the BAO is 
that we obtain the growth parameters at a higher redshift, which is complementary to the 
SNe results. In the future, if systematics are controlled, and probes like JDEM are able to 
measure both SNe and BAO data, we should be able to obtain independent estimates of 
the growth parameters at both very low and very high redshifts from this method. 
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4. Conclusions 

In this paper we have proposed a method for extracting growth parameters for dark 
energy models (within the spatially flat FRW universe) from observations that map the 
background universe, such as measures of the luminosity distance or the angular diameter 
distance. The method is model independent and unbiased. For current data, the growth 
factor g(z) may be constrained to ~ 5% at z — 0.3, while the growth rate f(z) is constrained 
to ~ 6%. For future JDEM SNe data, we will be able to put constraints of the order of a 
few percent on the growth parameters, e.g. 2% on the growth factor and 3% on the growth 
rate at a redshift of 0.3 , and 4% on the growth factor and 8% on the growth rate at a 
redshift of unity. In conjunction with the likelihood parameter estimation method, this 
method acts as an important consistency check on the accuracy of the priors on for 
SNe. With future probes like JDEM and BOSS taken in conjunction, it will lead to an 
unbiased estimation of the growth parameters upto a redshift of z = 2.5. 

It is well known that, in GR and for most DE models, the exp ansion history com pletely 



determines the linearized growth rate of density perturbations (jStarobinsky 1998 



Sahni fc Starobinsky 20061 ) (the exact conditions for this are formulated at the beginning 
of Sec. II). Consequently, a comparison of the density contrast reconstructed from the 
expansion history would provide one more important consistency check for a large variety 
of DE models including the cosmological constant and quintessence. On the other hand, as 
explained in more detail at the beginning of Sec. II, any departure of the observed density 
contrast from that reconstructed using standard candles and rulers would almost certainly 
indicate that either there is an exchange of energy between DE and DM (so that the effective 
energy-momentum tensor of DE is not on its own covariantly conserved), or that cosmic 
acceleration is a consequence of modified, non-Einsteinian gravity. In modified gravity 
theories, such as Braneworld models, scalar-tensor and f(R) gravity, etc. , the linearized 



perturbation equation for 5 does not follow the Newtonian form, Eq (j3J) (ILue et al. 2004 
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Kovama fc Maartens 20061; IShtanov et al. 20071; iHu fc Sawicki 20071: iBoisseau et al. 2000 



Jain fc Zhang 2007 



Song fc Kovama 2008 



Zhao et al. 2008 



Song fc Pore 20081 ). Hence 



the density contrast reconstructed using observations of standard candles/rulers via Eq f[T2T) 
and the density contrast determined directly from observations of large scale structure, 



say, by weak 



(IGuzzo 2008 



lensing, ga 



Song 2008 



axy redshift distortions or cluster abundances at different z 

1 I 1 

Vikhlinin 2008 : bimatti et al. 20081 ). are likely to differ. 



In Sec. 13.41 we show that one can obtain a strong signature of modified gravity by 
comparing results of this reconstruction method with future observations of galaxy redshift 
distortions using the DGP model as a toy example of modified gravity, where the growth 
factor g(z) is scale-independent on small scales. However, as discussed in Sec. II, g(z) 
often becomes scale-dependent both in modified gravity and in the case of direct DE-DM 
interaction (or their unification). Therefore, for further discrimination of DE models 
alternative to quintessence and the cosmological constant, measurement of 5 at different 
comoving scales is required to determine if g(z) is scale-dependent or not. 

Future surveys such as JDEM are expected to deliver high quality data for both 
supernovae and weak lensing. Using such surveys it would then be possible to compare 
the reconstructed density contrast from standard candles (SNe) with the density contrast 
observed from gravitational clustering (lensing). Therefore, we hope that the techniques 
developed in this paper, combined with future observations, will help unravel the nature of 
that most enigmatic quantity - dark energy. 
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Fig. 1. — Reconstructed linear growth factor g{z) for different datasets. The top panels show the 
results for Model 1 (ACDM) using Union-like (set A, left panel) and JDEM-like (set B, right panel) 
SNe datasets, while the bottom panels show results for Model 2 (variable w, eq (fTKj) ) using set A 
(left panel) and set B (right panel). In each figure, the black dotted line represents the true model, 
while the green dashed line represents the other model. The red solid lines show the la error bars 
for the integral reconstruction using eq (j 12 [> . 
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Fig. 2. — Reconstructed growth rate f{z) for different datasets. The top panels show the results 
for Model 1 (ACDM) using Union-like (set A, left panel) and JDEM-like (set B, right panel) SNe 
datasets, while the bottom panels show results for Model 2 (variable w, eq (|18p ) using set A (left 
panel) and set B (right panel). In each figure, the thick black dotted line represents the true model, 
while the green dashed line represents the other model. The red solid lines show the la error bars 
for the integral reconstruction using eq (|12p T he blue vertical line s in the right panel show the 



expected observational constraints from Euclid 



Cimatti et al. 2008 ). 
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Fig. 3. — Reconstructed growth parameters for different datasets using the smoothing scheme 
Eq (|16p on the integral reconstruction method, eq f)12[) . The top panels show the results for Model 
1 (ACDM) for the growth factor g(z) (left panel) and the growth rate f(z) (right panel). The 
bottom panels show the results for Model 2 (variable w, eq (|18|) ) for g(z) (left panel) and f(z) 
(right panel). In each figure, the black dotted line represents the true model, while the pink dashed 
line represents the other model. The green dashed shaded area represents the la errors for the 
integral reconstruction of set A (Union- like) , while the green hatched shaded area represents the 
reconstruction for set B (JDEM-like). The blue vertical lin es in the right panel show the expected 
observational constraints from Euclid (jCimatti et al. 20081 ). 
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Fig. 4. — Reconstructed growth rate f(z) for model 1 (ACDM) using various redshift distributions. 
We use (a) set A (Union-like) with number of SNe doubled at low and high redshifts (left panel) 
(b) set A with number of supernova doubled for mid-range SNe (center panel) and (c) JDEM-like 
(set B) redshift distribution with Union- like (set A) errors (right panel). In each panel, the red 
solid lines depict la error bars on set A, while the green dashed lines show the la error bars on set 
A modified according to (a), (b), (c). The black dotted line represents the true model. 
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Fig. 5. — Reconstructed growth rate f{z) for model 1 (ACDM) using various error distributions. 
We use (a) set A (Union-like) with errors halved at low and high redshifts (left panel) (b) set A 
with errors halved for mid-range SNe (center panel) and (c) set A with JDEM-like errors for each 
SNe (right panel). In each panel, the red solid lines depict la error bars on set A, while the green 
dashed lines show the la error bars on set A modified according to (a), (b), (c). The black dotted 
line represents the true model. 
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Fig. 6. — Reconstructed growth rate f(z) for model 1 (left panel) and model 2 (right panel) using 
set B (JDEM-like) with different reconstruction methods. The red solid lines show the la limits 
for reconstructed f(z) using the integral reconstruction method, eq (|12p . while the green hatched 
region shows the la limits for f{z) using w parameterization, eqs (|19|) . (|21|) . The black dotted line 
represents the true model. 
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Fig. 7. — Reconstructed growth rate f(z) for model 1 (left panel) and model 2 (right panel) for 
set B (JDEM-like) with different reconstruction methods, using f2om = ^Om(t r ue) + 0.03. The 
red solid lines show the la limits for reconstructed f(z) using the integral reconstruction method, 
eq (|12p . while the green hatched region shows the la limits for f(z) using w parameterization, 
eqs (I19D . (j21j) . The black dotted line represents the true model. Note that the results for the two 
reconstructions lie on opposite sides of the true value of f(z). 
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Fig. 8. — Reconstructed growth rate f{z) for JDEM-like dataset using a modified gravity model 
(DGP, eq (|23p ). The thick black dotted line represents the result expected from just the expansion 
history (eq (fT9j) using eq (|23|) ). while the green dashed line represents the result expected from 
gravitational clustering (eq (|24p ). The red solid lines show the la error bars for the integral 
reconstructio n using eq (I12I). Th e blue vertical lines show the expected observational constraints 



from Euclid (jCimatti et al. 2008l ). The discrepancy between the two would act as a signal for 



modified gravity. 
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Fig. 9. — Reconstructed growth rate f(z) for current Union set of SNe data, using Oom = 0.26 ± 
0.03. The red solid lines show the la limits for reconstructed growth parameter using the integral 
reconstruction method, eq (|12|) . while the green dashed shaded area shows the la limits for the 
parameter using smoothing scheme, eq (|16p . for the integral reconstruction methods. The black 
dotted line shows f(z) for LCDM, the green dashed line shows f(z) for Model 2 (variable w, 



eq (fT8l)l. The th ree vertica" 



blue lines show the current measureme nts of f(z) from 2dFGRS 



(jVerde et al. 20021 ). 2SLAQ (jRoss et al 20071 ) and VVDS (jGuzzo 2008h . 
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Table 1: Reconstructed linear growth factor g and growth rate / using different datasets for Model 1 



Datasets 


z 


9{z) 


^/smooth (^) 


<?exact(-^) 


m 


/smooth (^) 


/exact (•£•) 




0.3 


1.11 ±0.04 


1.12 ±0.02 


1.12 


0.65 ±0.04 


0.63 ±0.03 


0.64 




0.6 


1.18 ±0.05 


1.18 ±0.03 


1.19 


0.77 ±0.06 


0.75 ±0.04 


0.76 


A (Union SNe) 


1.0 


1.27 ±0.06 


1.26 ±0.05 


1.25 


0.82 ±0.09 


0.83 ±0.06 


0.85 




1.5 


1.24 ±0.15 


1.26 ±0.09 


1.28 


0.97 ±0.21 


0.94 ±0.10 


0.92 




0.3 


1.13 ±0.02 


1.12 ±0.01 


1.12 


0.64 ±0.02 


0.63 ±0.01 


0.64 




0.6 


1.21 ±0.03 


1.20 ±0.02 


1.19 


0.75 ±0.04 


0.76 ±0.02 


0.76 


B (JDEM SNe) 


1.0 


1.24 ±0.05 


1.23 ±0.03 


1.25 


0.86 ±0.07 


0.84 ±0.03 


0.85 




1.5 


1.25 ±0.08 


1.26 ±0.09 


1.28 


0.93 ±0.11 


0.92 ±0.06 


0.92 


C (BOSS BAO) 


2.5 


1.28 ±0.13 


1.29 ±0.11 


1.30 


1.01 ±0.09 


1.00 ±0.07 


0.97 



Table 2: Reconstructed linear growth factor g and growth rate / using different datasets for Model 2 



Datasets 


z 


9{z) 


^/smooth (^) 


<7exact 


m 


/smooth 


/exact 




0.3 


1.12 ±0.03 


1.12 ±0.01 


1.13 


0.60 ±0.04 


0.61 ±0.03 


0.61 




0.6 


1.20 ±0.05 


1.20 ±0.03 


1.21 


0.70 ±0.05 


0.72 ±0.05 


0.72 


A (Union SNe) 


1.0 


1.26 ±0.07 


1.27 ±0.05 


1.28 


0.82 ±0.09 


0.81 ±0.07 


0.80 




1.5 


1.34 ±0.18 


1.33 ±0.10 


1.32 


0.89 ±0.20 


0.90 ±0.13 


0.87 




0.3 


1.13 ±0.02 


1.11 ±0.01 


1.13 


0.62 ±0.02 


0.61 ±0.01 


0.61 




0.6 


1.19 ±0.04 


1.20 ±0.02 


1.21 


0.71 ±0.03 


0.70 ±0.02 


0.72 


B (JDEM SNe) 


1.0 


1.27 ±0.05 


1.27 ±0.04 


1.28 


0.81 ±0.04 


0.80 ±0.03 


0.80 




1.5 


1.33 ±0.08 


1.31 ±0.07 


1.32 


0.88 ±0.08 


0.86 ±0.04 


0.87 


C (BOSS BAO) 


2.5 


1.37 ±0.13 


1.35 ±0.07 


1.34 


0.96 ±0.11 


0.94 ±0.06 


0.93 
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Table 3: Reconstructed linear growth factor g and growth rate / using current supernova data 



z 


sO) 


^smooth (^) 




/(*) 


/smooth (-2*) 


/acdm(z) 


0.3 


1.13 ±0.05 


1.13 ±0.05 


1.12 


0.62 ±0.06 


0.61 ±0.04 


0.64 


0.6 


1.21 ± 0.07 


1.20 ±0.06 


1.19 


0.79 ±0.08 


0.80 ±0.05 


0.76 


1.0 


1.29 ±0.10 


1.28 ±0.07 


1.25 


0.93 ±0.11 


0.88 ±0.08 


0.85 


1.5 


1.37 ±0.16 


1.35 ±0.12 


1.28 


1.05 ±0.24 


0.98 ±0.11 


0.92 



